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Abstract 

We study the off-equilibrium relaxational dynamics of the Amit-Roginsky 
3 field theory, for which the mode coupling approximation is exact. We show 
that complex phenomena such as aging and ergodicity breaking are present at 
low temperature, similarly to what is found in long range spin glasses. This 
is a generalization of mode coupling theory of the structural glass transition 
to off-equilibrium situations. 
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Recently there has been important progress in the understanding of the off- 
equilibrium dynamics of disordered systems. Although aging phenomena were dis- 
covered experimentally in spin glasses in 1983 IJ and subsequently studied by several 
groups in spin glasses [Q] and many other different disordered systems, only quite 
recently has a dynamical mean field theory (DMFT) for non-equilibrium phenom- 
ena, based on microscopic spin glass models, been put forward P, f§, ||. Among the 
most interesting feaatures which emerge from the theory are the violations at low 
temperature of time translation invariance (TTI) and of the fluctuation-dissipation 
theorem (FDT), even for very large times, implying breaking of ergodicity. This 
breaking of ergodicity can not be described simply as a separation of the phase 
space into different ergodic components within which the system thermalizes as 
in usual broken-symmetry phases. Rather, the distribution function continues to 
spread over larger regions of configuration space as time goes on. This spreading 
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gets slower and slower but never stops; thus the system never comes to equilibrium, 
and its properties depend on the time elapsed since the quench to low temperature. 

Far from being confined to the realm of disordered system aging phenomena have 
been observed in a number of different systems without disorder, in experiments |J 
and numerical simulations [^]. 

At present there is a great deal of theoretical research on spin-glass like phe- 
nomena in non-random systems || || 10 1. In || || a mechanism has been proposed 



in which an effective disorder can be induced in a priori uniform systems (althoug 
with "complicated" interactions among the variables). In this paper we focus on 
the dynamical mean field theory for glasses, usually called mode coupling theory 
(MCT), showing how it is possible to extend it in order to include non equilibrium 
phenomena. 

It has been noted |nj that the mean field dynamical equations of many disor- 
dered mean field models have, in the high temperature phase, a structure which is 
surprisingly similar to that of the equation of the mode coupling theory often used 
to describe the structural glass transition |L2|, |13| . The mode coupling equations 
have been derived from molecular or hydrodynamical models based on an approx- 
imation in which vertex correction are neglected in the perturbation expansion of 
the self-energy jl4H . Systematic use is made of the hypothesis that the system under 
study is at thermal equilibrium. On the other hand there is experimental evidence 
that real glasses are out of equilibrium [H| . Even in MCT the ideal glass transition 
is described as one into a non-ergodic state. Thus current mode coupling theory can 
at best give a good description of the physics at high enough temperature. 

This letter is a step towards the extension of the MCT below the glassy temper- 
ature, including the possibility of studying nonequilibrium phenomena. Technically, 
we construct a theory in which we do not assume TTI or the FDT ah initio in the 
derivation of the equations. As an example we have chosen to study a simple model 
without quenched disorder, where the mode coupling approximation is exact. This 
is the Amit-Roginsky </> 3 field theory |TJ|, originally proposed to study the critical 
properties of the 3-state Potts model. We will find dynamical equations identical to 



that of the spherical j>-spin glass model [16, 11, p[ for p = 3. At high temperature 



supposing TTI and FDT one gets usual MCT equations, similar to those first pro- 
posed by Leutheusser to study the glassy transition. At low temperature there is 
a phase where ergodicity is broken in a way now known from mean field spin glass 
dynamics. 

Amit and Roginsky (AR) consider an 0(3) symmetric 3 field theory where the 
basic fields transform under rotation according to representations of high angular 
momentum quantum number I. Recently the same method has been exploited in 



two interesting papers |17fl , to find soluble limits of nonlinear stocastic equations. 
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The method is applicable any time the nonlinearity is represented by a bilinear term. 
In this letter we discuss for simplicity the dynamics of the original AR model. In 



a forthcoming paper |L8| we will discuss the case of the fluctuating compressible 



hydrodynamics of Das and Mazenko |T3|], which, although formally much more com- 
plicated than the present model, gives rise to equations of with the same structure. 
The AR model is defined by the 0(3) invariant Lagrangian 
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with iV = 21 + 1, and ip m transforms according to the irreducible N— dimensional 
representation of 0(3). In the following we will denote by Ci[ip,ip*} the cubic part 

of the lagrangian. The coefficients f ' | are the Wigner 3?'-symbols. 

y mi m 2 m 3 J 

To have a non vanishing interaction term / must be even. Amit and Roginsky have 
shown that in the limit / — > oo vertex corrections in the perturbation expantion 
can be neglected. Thus, defining the equilibrium corrlation function G(x, y) = 
Yjm{ % l J m{,x) 1 Pm(y)) 'GMs, the perturbative series can be resummed, and the Dyson 
equation of the model in real space is simply: 

G-\x,y) = G 1 (x,y)- 9 j[G(x,y)} 2 , (2) 

where Gq (x, y) = 5(x — y)(fi 2 — V 2 ). Similarly one finds that the free energy as a 
functional of G can be written 

F = — TrLog(G) + J dxdyG^(x, y)G(y, x) - ^ J dxdyG(x, yf (3) 

and (H) is its associated variational equation. It is assumed in deriving ([|) that 
(Wn(aO) = 0- 

In order to avoid problems connected with the unnormalizability of the Gibbs 
measure for a 3 theory, we will consider a "spherical" variation of the AR model, 
imposing the constraint Xim'/'mfilCi 1 ) = Nq. In this situation, /i 2 can be viewed 
as a Lagrange multiplier which enforces the constraint. Eq. (0) for x = y should in 
this case be read as an equation for /z 2 . 

We study the relaxational dynamics of this model, restricting our attention to 
the zero-dimensional case (we drop the x dependence and the gradient term in ([I])). 
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The extension of our results to finite dimension is immediate. We consider the 
Langevin equation 

-1> m {t) = -ii 2 (t)if> m (t) - + rj m (t) (4) 

where rj m is a white noise with correlation function (Vm(t)Vn(. s )) = 2T5 mn 8(t — s), 
and T is the temperature. We will see a posteriori that it is possible to choose [i as 
a function of t to impose the spherical constraint at all times. In order to generate 
the perturbation expansion of eq.([|), we consider the Martin-Siggia-Rose generating 
functional 

Z[h] = J Vipvr (j[ 6 (j t ^ m (t) + - Vm(t)j ) exp (| dth* m (t)^ m (t)) JbP] 

(5) 

J\}\)\ is the functional determinant of d/dt + d 2 C t / difjmdip'^ which, with the Ito 
convention, is equal to 1, and we denoted by C t the time dependent "Lagrangian" 
C t = fi 2 (t) Ysm i : mi J m+£i- We now follow a standard procedure |19| to prove that the 



structure of the diagrammatic expansion for dynamics is the same as that for statics. 
We introduce an auxiliary field iip m (t) } to get Z = J UipVip*UipVi/j* exp / dtitp^ 

[ip m — iipm + (dC/dtp^)}^ . For simplicity we have set h m = 0. This will not affect 
our subsequent manipulations. We then introduce an even grassman variable 8, and 
the "superfield" (fi m (t) = ip m (t) + 0itp m {t). In terms of these, Z reads, 

VWexp (-fat J d#{CM)[(i-^)^-T^^^ 

(6) 

The time integration here should be intended to start from an initial time t — 0, 
where the system is in a random initial condition. This corresponds to the experi- 
mental situation of a sudden quench of the system from very high temperature at 
the initial time. 

The interaction term in (|6|) is represented by the same function Cj as in statics, 
but its arguments are now the superfields (f> and (f>* rather than just of ip and ip*. Thus 
the structure of the dynamical perturbation expansion in terms of the superfield is 
formally identical to that of statics in terms of the field. We define at this point the 
correlator 

G(t,6;t',6>) = (0 m M)C(^')> 

+ e'^ut^Ut')} + oo'^UtW m (t')}. (7) 
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We see that G codes for the correlation function C(t,t') = (■ip m (t)ip^ n (t')) and the 
response function r(t,t') = (ipm(t)iipm(t')) = (5ip m (t) / Srj m (t')) , while the conserva- 
tion of probability implies {iipmit^ip^it')) = 0. Causality and the Ito convention 
imply r(t, t') = for t' > t and r(t + , t) = 1. The zeroth order correlator in dynamics 
is defined by Gq 1 = (1 — J7 — TO-jL + ^ 2 (t). In terms of G and G defined in 
this way, the Dyson equation of the theory has the same form as (0), replacinging 
x and y respectively by (t, 6) and (t',6'). The appearence of C and r only in the 
combination given by G is a consequence of the gradient character of the Langevin 
equation. 

Inverting the equation (H), multipling it by G, and disentangling the superfield 
notation, we get coupled equations for the correlation and the response function; for 
t > t' we have: 



-fi 2 (t)C{t,t') + ]-g 2 \ f ds2C{t,s)r{t,s)C{t',s)+ f dsC 2 {t, s)r(t' , s) 
2 [jo jo 

-V 2 (t)r(t 7 t') + lg 2 [ t ds2C(t,s)r(t,s)r(s,t') ( 
z Jo 



dC(t, t') 

dt 
dr(t,t') 
dt 

while for arbitrary fi 2 (t), C(t,t) satisfies: 



As mentioned above, we can choose to satisfy J2mi J m' l Prn = Nq at all times. We 
recognize in (ROD well studied equations in mean field spin glass theory: they are 



identical to those found in the dynamics of the so called spherical p-spin model [11] 
for p = 3. This is defined by the hamiltonian 

H — — 2_j JijkSiSjSk (10) 

i<j<k 

where the "spins" 5"s are real veriables satisfing the constraint J2{Li S 2 = Nq, and 
the couplings Jijk are independent Gaussian variables, symmetric under interchange 
of any pair of indices, with zero mean and variance Jfj k = 3g 2 /N 2 . The equivalence 
for large iV of the 0(3) symmetric vertex with a random one in dynamics had already 
been noted in fTTH . Note that eq. (|,pj) are causal first order integro-differential 
equations which admit a unique solution for finite times. 

In a seminal paper || Cugliandolo and Kurchan have found an asymptotic so- 
lution to the equation of the p-spin model, showing the existence of a low temper- 
ature phase in which TTI and FDT are violated even in the infinite time limit. 
The explicit numerical integration of equations with the same structure obtained 



5 



for a related model |IJ] support the asymptotic analysis of 0, which we review 
briefly here in the present context. At high temperature, the asymptotic dynamics 
is completely characterized by the two functions C as {r) = limt->ooC(t + r,t) and 
r asij) = H'm>t-+oo r (t + r,t) related by FDT, and such that lim,--^ C as {r) = 0. 
The long time limit of \x is given by /i 2 = (g 2 /2T)q 2 + T/q. This is just the 
zero- dimensional limit of eq. (Q) for the statics. Accordingly the "energy" of 
the system, E = (Cj(t)) = — (g 2 /3T) J * dsC 2 (t, s)r(t, s) tends to the static value 
E static = — 9 2 q 3 /(3T). The set of equations (§,0) reduces to a single one e.g. for the 
response function, almost identical to the one originally proposed by Lautheusser 



Below the critical temperature T c = g{q/2)^/ 2 \ such a regime still exists for 
t, t' — > oo with r/t and r/t' — > 0, but with lim^oo C as (r) = q ^ 0. In addition, 
there is an aging regime, in which TTI and FDT are violated. In this regime the 
limit t,t' — > oo is taken fixing the ratio h(t')/h(t) = e~ a (cr > 0).[] Here one has 

C(t,t') = C{a) 

Continuity between the homogeneous and aging regimes, implies C(0) = q. The 
value of q is specified at low temperature by the largest root of the equation 

T 2 1 

q=-,T~~ v? ( 12 ) 



g 2 (q - q) 2 



while jj 2 and u are given by 



f = ^9\f ~ ? 2 ) + j^—y u = L-l (13) 

2T (q-q) q 

It is worth noticing that the former equations can be found in the replica formalism 
for the £>-spin glass by imposing the condition of marginal stability. The energy takes 
contribution both from the asymptotic and the aging regins of times, and tends to the 
limit Eoyn = —g 2 (q 3 — (1 ~ M )<? 3 )/3T, different from its static counterpart. Figure 
1 shows the results of numerical integration of eqs. (H) and @ at a temperature 
where aging occurs. 

We have seen that the dynamics of the spherical AR model quite unexpectly 
displays glassy behaviour, with a sharp phase transition from ergodic to non-ergodic 



1 h(t) is an arbitrary increasing function which the present theory is not able to predict. On a 
numerical basis it was claimed in Q that h(t) — t, leading to the suggestive result C(t, t) — C(t'/t). 
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behaviour. For a correct picture of the model at low temperature, a set of coupled 
equations for the correlation and response functions have to be solved. It has been 
often discussed in the literature whether glasses undergo a real phase transition, 
or whether a better picture is one of a progressive freezing |TJ], MC equations 
which fit better with this scenario have been proposed [14, 13]. Even though in this 
case one can not expect aging behaviour for infinite times no matter how the limit 
is taken, the relaxation to equilibrium can be extremely slow. On an experimental 
time scale aging phenomena (or eventually "interrupted aging") can be present, and 
the use of the off-equilibrium equations can be more appropriate than that of the 
asymptotic one. 

A natural generalization of the MC eq. (H), and consequently of is obtained 
by replacing the non linear term by some generic "self-energy" function H(G(x,y)). 
This can be introduced on a purely phenomenological ground, as it has been done 
by Gotze in the equilibrium case ||20|| , or derived from some microscopic theory. 
Examples of equations of this form have been analyzed extensively in [5] in the 
context of mean field disordered systems. Depending on the form of S, the long time 
aging regime can assume different forms. In the model discussed here, we have found 
the simplest among the scenarios proposed in ||, ||. Further work is necessary to 
determine the most appropriate form of S to interpret real glass experiments. Our 
results however, open the door for the systematic study, starting from a mean-field 
level, of aging effects in real glasses. 



It is a pleasure to thank C. de Dominicis, M. Mezard, G. Parisi for interesting 
discussions and suggestions. 
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Figure 1: 



Figure caption 

Figure I. A sketch of the correlation and the response functions for g = q = 1, 
T = .2 < T c . (a). Normalized correlation function C(t w + r,t w )/q and integrated 
response function m(t w + r,t w ) = (T/q) Jq w r(t w + r, s)ds as functions of r, for 
different waiting times t w . If TTI and the FDT held, these functions would be equal 
up to a constant shift and independent of t w . (b). m plotted against C for different 
waiting times. For large values of C (short times), the FDT holds, and the curves 
approach a straight line of unit slope which intersects the C axis at q. The deviation 
from this line at smaller C (long times) indicates off-equilibrium behaviour. In the 
region C < q the behaviour predicted by the asymptotic solution would be a straight 
line with slope u, starting from the origin. The times displayed in the figure are 
very short compared to the ones needed to see the asymptotic behaviour. 
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